Ground states of hard-core bosons in one dimensional periodic potentials 
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With Girardeau's Fermi-Bose mapping, we find the exact ground states of hard-core bosons 
residing in a one dimensional periodic potential. The analysis of these ground states shows that 
when the number of bosons A*' is commensurate with the number of wells M in the periodic potential, 
the boson system is a Mott insulator whose energy gap, however, is given by the single-particle band 
gap of the periodic potential; when N is not commensurate with M, the system is a metal (not 
a superfluid). In fact, we argue that there may be no superfluid phase for any one-dimensional 
boson system in terms of Landau's criterion of superfluidity. The Kronig-Penney potential is used 
to illustrate our results. 



I. INTRODUCTION 

A Tonks- Girardeau (TG) gas is a one dimensional 
strongly correlated system consisting of bosons with 
hard-core interaction. This model of a one dimensional 
gas was first studied by Tonks 1], who treated it as a 
classical gas. Later studies were nevertheless focused 
mostly on its properties as a quantum boson gas[2,B- As 
the milestone development in this model, Girardeau^, Q 
found that this boson gas can be mapped into a spinless 
fermion system and obtained all of its eigenstates. 

This idealized TG gas had only been a curiosity of 
theorists until the realization of Bose-Einstein conden- 
sation with dilute atomic gases in 1995 [1, 0, It was 
soon noticed that the TG gas could be realized exper- 
imentally with a Bose-Einstein condensate trapped in 
tightly-confined waveguides Q. This work renewed inter- 
ests in the TG gas and led to a series of further theoret- 
ical studies^, [nl, [H [H, [ii, [IE [H, [13 ■ In particular, 
Lapeyre et al. studied the momentum distribution of a 
harmonically trapped TG gas [11] and Forrester [l^ et al. 
obtained highly accurate results on the TG gas's long- 
range off-diagonal behavior [Tsl [ii]. Finally in 2004, the 
TG gas was first realized with a Bose-Einstein conden- 
sate trapped in two perpendicular optical latticespol.l2l|. 
An overview of studies in this interesting subject can be 
found in a recent review by Yukalov and Girardeau [2^. 

In this article we study the TG gas in a periodic poten- 
tial. We focus on the case that the hard-core interaction 
between bosons is a repulsive infinite contact interaction. 
By applying Girardeau's Fermi-Bose mapping^, Q, this 
Boson gas can be mapped into a spinless free Fermi sys- 
tem. As a result, we can obtain all the eigenstates and 
eigen-energies of this TG gas. We are particularly inter- 
ested in its ground properties. Our analysis of the ground 
states shows that such a TG gas has two quantum phases, 
Mott insulator and metal. When the number of bosons 
N in the gas is commensurate with the number of wells 
M in the periodic potential, the boson system is a Mott 
insulator; however, its energy gap is given by the single- 
particle band gap of the periodic potential. When N is 
not commensurate with M, the system is a metal. We 
emphasize that this boson metal is not a superfluid as 
it can not support any superflow with finite critical ve- 



locity. In fact, in this sense there may be no superfluid 
phase for any one-dimensional boson systems including 
soft-core bosons. 

The Kronig-Penney potential is used to illustrate our 
results. Various properties of the ground state are com- 
puted, such as pair distribution function, single-particle 
density matrix, and momentum distribution. Before we 
proceed with full discussion, we note that the hard- 
core bosons in a lattice have been actively studied un- 
der single-band approximation with the Jordan- Wigner 
transformation 23, 24, 25, 2^- It is also studied as a 
limiting case in Ref . |27l| 



II. TG GAS IN PERIODIC POTENTIAL 

Normally, the boson in a TG gas has an "impenetra- 
ble" hard core characterized by a radius of a. In this 
study, we focus on the case that the hard core is a point 
with no radius, that is, the inter-particle interaction is 
given by 



U{x) 



oo 




x = 0, 



(1) 



Such an interaction is equivalent to a constraint on the 
wave function ip{xi, • • • , xn)- 



= 



if 



Xj, 1 ^ i < j ^ N. 



(2) 



This means that this TG gas can be viewed as a group 
of "free" bosons governed by the following Hamiltonian 



H 



2m dx-j 



(3) 



while its wavefunction is subject to the constraint in 
Eq.(I2]). Based on the observation that the constraint in 
Eq.(I2]) is automatically satisfied by any wavefunction of 
a Fermi system due to its antisymmetry, Girardeau 
found a natural mapping that allows one to construct 
wavefunctions for a TG gas from wavefunctions of a free 
Fermi system. The mapping is 



ip^{xi,X2,-- ■ ,xn) = Aip^{xi, 



X2, 



, Xn) (4) 
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where i/;^ is the Bose wavefunction and tp^ the Fermi 
wavefunction. A is caUed unit antisymmetry function 
which is defined as 



A{xi,X2, ■ ■ ■ ,xn) = ]^sgn(a;i - Xj) , 



(5) 



i>j 



where sgn is the sign function. This Fermi-Bose mapping 
is one-to-one and, therefore, the energy spectrum of a 
TG gas of infinite contact potential is the same as a free 
spinless Fermi system0,Q. 

In this article, we study the case where the external 
potential V{x) is periodic, V{x) = V{x + d), with d be- 
ing the period. Due to the Fermi-Bose mapping we 
consider first the system of N free spinless fermions re- 
siding in this periodic potential. According to the basic 
knowledge of solid state physics [1^, when the number of 
fermions N in the gas is commensurate with the number 
of wells M in the periodic potential, the fermions can fill 
up exactly N/M Bloch bands and thus the Fermi system 
is an insulator with its energy gap determined by the pe- 
riodic potential. When N is not commensurate with M, 
the fermions can fill up the bands only partially and the 
Fermi system is a metal. 

Because of the Fermi-Bose mapping (j?]), the TG gas 
has the same energy spectrum as this free Fermi system. 
This allows us to conclude immediately that when N is 
commensurate with M, the TG gas is an insulator. How- 
ever, we emphasize that this insulator is a Mott insulator 
unlike its mapping target, the free Fermi system, where 
the insulator is a band insulator. This leads to a quite 
peculiar situation, the energy gap of a Mott insulator is 
completely dictated by the periodic potential and is given 
by the single-particle band gap. For the other situation 
where N is not commensurate with M, the TG gas is a 
metal. 

We emphasize that this Bose metal is a real metal 
not a superfiuid according to Landau's criterion of 
superfluidity d^. Landau's criterion is that the system is 
a superfluid if it has phonons as its only low energy exci- 
tations. Plotted in Fig[T]is the excitation spectrum of a 
TG gas in a periodic potential for the non-commensurate 
case, which is marked out by the shade. It is clear from 
this figure that there are two kinds of low energy excita- 
tions, one at fc = and the other at fc = fcmax- The low 
energy excitation at fc = may be barely called phonon 
since the phonon relation e — ck exists only at the limit 
fc — > and the excitation spectrum immediately spreads 
out at fc 7^ 0. The most destructive is the low energy ex- 
citations at fcmax, which are not phonons. The presence 
of these low excitations make it vanishing the minimum 
slope of e{k)/k, which gives the critical velocity of pos- 
sible superfluidity. Therefore, we can conclude that this 
TG gas can not support a superflow with finite critical 
velocity and it is not a superfluid. 

What is more interesting is that the TG gas is not an 
isolated case. According to Ref.[3^, a boson gas with 
finite repulsive ^-function interaction has a very similar 
excitation spectrum as our TG gas. This means that such 



a gas is no superfiuid, either. This prompts us to spec- 
ulate any one-dimensional boson system is not a super- 
fiuid. We are aware that there are different definitions of 
superfiuiditvfsH . In particular, according to Leggettfs^, 
our TG gas (in his version it is spinless fermion) is a su- 
perfluid since it has the so-called non-classical rotational 
inertia (NCRI). Because of this, we state again that we 
are talking about superfluid in the sense of Landau's cri- 
terion. 
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FIG. 1: Excitation spectrum of a TG gas in a periodic poten- 
tial for the non-commensurate case. The shade makes out all 
the possible excitations, fcmax is the largest possible value of 
fc, similar to the Fermi wavevector in the fermion case. 

We now construct explicitly the exact eigenfunctions 
of this TG gas. The single-particle eigenstate (pn,k{x) is 
a Bloch state that satisfies the Schrodinger equation 

^"2to 9a?2 ^i^)\fn,k{x) = E„^k(PnM^) ■ (6) 

The notation n is for band index and k is the Bloch wave 
vector in the first Brillouin zone. For convenience, we 
use a = {n, k} to denote both the band index and the 
Bloch wave vector. Then according to Eq.Q, the Bose 
eigenfunction is given by the Slater determinant 



A 



,X2r ■ ■ ,xn) = 

Vaiix2) Va2{x2) ■■■ Va„(a;2) 
^ai{xN) ^uAxn) ■■■ 'PaNixN) 



(7) 



As one can check, this symmetric wavefunction satisfies 
both the Schrodinger equation with the Hamiltonian Q 
and the constraint ([2]). The corresponding eigen-energy 
is 



N 



(8) 
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For the ground state, the N eigenfunctions in the Slater 
determinant are for the N lowest eigenstates. For conve- 
nience, we shall assume that both N and M are odd in 
the following discussion. 



III. TG GAS IN THE KRONIG-PENNEY 
POTENTIAL 
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FIG. 2: Width of the energy gap between the first and the 
second Bloch band as a function of the strength V of periodic 
potential. This is also the energy gap of the Mott insulator 
when the number of bosons in the TG gas is equal to the 
number of wells in the periodic potential. 



We now use the Kronig-Penney (KP) potential to illus- 
trate our general results of the TG gas in the last section. 
The Kronig-Penney potential is given by 



(9) 



where 7 is the strength of the J-potential. We follow the 
standard procedure in solid state physics [1^ by placing 
the system in a box whose size is L = Md and impos- 
ing the usual periodic boundary condition. If we use d, 
the period of the potential, as the unit of distance and 
2m(P /h^ as the unit of energy, the Schrodinger equation 
can be written as 



d 



2 M-l 



[^g^ + Vj2Si^-j)Mx)^E^^^{x) (10) 



3=0 



where V = 2m'yd/h?. The eigenfunction ipa is given by 



fix) = A„e'P°" + B„e-'P°" X e [0, 1) 
M^) = { e'^^'fix-s) xe[s,.s + l) (11) 

(s = l,2,.-. ,M-1) 



with the coefficients determined by 



27r/ , M -1 

V . 

cos Pa H Smpa — COS K 

2pa 



Br, 



y^e-'P^/^Pa sm[{pa + k)/2] (12) 
^/[{2pl + V) sinpa - PaV cos Pa]M sin Pa 

V2e'P-^^Pa sin[(p„ - fc)/2] 
y/[{2pl + V) sinpa - PaV cospa]M sinpa 



The corresponding eigen-energy is Ea = p^ ■ 

As we have demonstrated in the last section, when a 
TG gas of exactly M bosons is placed in the KP potential, 
the system is a Mott insulator. The energy gap between 
the first and the second Bloch band is plotted in Fig.l. 
This figure shows that this gap initially increases with 
the lattice strength V and eventually saturates at the 
energy difference between the first excited state and the 
ground state in a square well potential. 



A. PAIR DISTRIBUTION FUNCTION 

The pair distribution function, normalized to N(N 
1), is defined as 



D{xi,X2) = NiN -1) / \i;^{xi,--- , XN)\^dx3 ■ ■ ■ dx 



N 



N-1 



(13) 



It is the joint probability of finding simultaneously one 
atom at xi and another atom at X2. When xi — X2, the 
pair distribution function vanishes by antisymmetry, re- 
flecting physically the impenetrable hard-core interaction 
between boson particles. 

Figure [3] shows gray-scale plots of the pair distribu- 
tion function D{xi, X2) for different numbers of particles 
of M^9. Some qualitative features of the pair distribu- 
tion function are observed. First, the figure shows that 
D{xi, X2) vanishes at contact xi = X2, reflecting the im- 
penetrability between the particles. Second, the peri- 
odicity due to the Kronig-Penney potential is apparent. 
Third, with the increase of the particle number, the black 
diagonal stripe becomes thinner. This indicates that the 
averaged distance between particles decreases. 



B. REDUCED SINGLE-PARTICLE DENSITY 
MATRIX 

The reduced single-particle density matrix is given by 
iv/<(x, 

^pQ {x' ,X2, ■ ■ ■ ,XN)dx2 ■ ■ ■ dXN- (14) 
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FIG. 3; Gray-scale plots of the pair distribution function 
D{xi,X2). V=l. (a) N^2, M=9; (b) iV=3, Af=9; (c) iV=6, 
M=9; (d) iV=9, M=9. 
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FIG. 4: Gray-scale plots of the reduced density matrix p(x, x') 
of the TG gas. V=l. (a) 7V=1, M=7; (b) 7V=3, M=7; (c) 
iV=5, M=7;(d) iV=7, M=7. 



The existence of the ofF-diagonal long-range order in this 
matrix indicates the onset of Bose-Einstein condensation 
[33I [131 • Its diagonal term p{x) = p{x,x' = x) is the 
single-particle density and is normalized to N, 



p{x,x)dx^N. (15) 
The multidimensional integral in Eq. p^ is evaluated 



numerically by Monte Carlo integration. The results are 
shown in Fig[4l The relative darker area of the diagonal 
in the figure is due to (5-potentials, which tend to repulse 
the particles away. The particles like to stay between 
(5-potentials, i.e., the brighter diagonal area. 
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FIG. 5: Off-diagonal elements p(x, M 
density matrix. V = I. (a) iV=l, M=7; 
iV=5, Af=7; (d) iV=7, M=7. 



- x) of the reduced 
(b) iV=3, M=7; (c) 



The gray scale in Figd] indicates that the ofF- 
diagonal elements decrease as A'' increases. To see this 
more clearly, we have plotted the off-diagonal elements 
p{x, M — x) of the reduced density matrix in Fig. [5] It is 
clear from the figure that as N increases and the interac- 
tion in the system gets stronger, the off-diagonal elements 
decrease relative to the diagonal. This indicates that the 
interaction suppresses the off-diagonal terms. It has been 
proved that there is no real Bose-Einstein condensate in 
the uniform Tonks- Girardeau gas p^. [l^ [Toj. 



C. MOMENTUM DISTRIBUTION 

The Bose-Einstein condensation of a Bose system is of- 
ten characterized by a macroscopic number of bosons oc- 
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FIG. 6: Normalized momentum distribution function Q{k). 
1/ = 1. M = 7. (a-d) Free Fermi gas for Af=l,3,5,7 (from top 
to bottom); (e-h) TG gas for A'^=l,3,5,7 (from top to bottom). 



cupying the zero momentum state. It is therefore mean- 
ingful to look at the momentum distribution of our TG 
gas in the ground state. The momentum distribution 
g{k)js related to the reduced density matrix and given 



g{k) 



1 



2ttN 



+ 00 



dx 



+ 00 



dx'p{x,x')e''^^'-'-'''\ (16) 



which is normalized to one. As noted by Girardeau [J, 
the Fermi-Bose mapping in Eq.Q becomes simply 
— for the ground state. This shows that even 

though the density distribution is the same for the 



TG gas and its mapping counterpart, free Fermi gas, their 
momentum distributions are different. It is therefore also 
interesting to compare the momentum distributions for 
these two different systems. 

In FiglGl we have plotted the momentum distributions 
for both the TG gas and the free Fermi gas in a peri- 
odic potential. The number of wells are M = 7 and the 
strength of the periodic potential \s V = 1. The left 
column (a-d) of Figl6] is for the Fermi gas and the right 
column (e-h) is for the TG gas. Since higher order Bragg 
peaks, e.g., the ones at ±27r, are very small, the plots in 
FigE] focus on the central Bragg peak, which consists of 
discrete small peaks due to the finite size of the systems. 

We first look at the right column of FiglS] For iV = 1, 
which can be regarded as the free boson case, we see that 
the Bragg peak is completely located at A: = 0. As the 
number of bosons increases, the interaction in the system 
gets stronger. This causes the Bragg peak to spread out 
as seen in the right column of Fig[6l When N = M, 
the peak has a width of almost a whole Brillouin zone. 
It is interesting to compare the TG gas with the free 
Fermi gas. As shown in the left column of FiglHl the 
Bragg peaks are much "fatter" than the TG gas (except 
iV = 1). This demonstrates that even a free Fermi gas 
has a broader momentum distribution than a strongest 
interacting boson gas. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have studied the ground-state prop- 
erties of a Tonks-Girardeau gas in a periodic potential 
with the Girardeau's Fermi-Bose mapping. We found 
that such a TG gas is a metal when the number of parti- 
cle N is not commensurate with the number of wells M 
in the periodic potential; it is a Mott insulator when N is 
commensurate with M. What is more interesting is the 
energy gap in the Mott insulator is given by the single- 
particle band gap of the periodic potential. To illustrate 
our results, we have employed the Kronig-Penney poten- 
tial to compute various properties of the ground state, 
such as pair distribution function, single-particle density 
matrix, and momentum distribution. 
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